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ABSTRACT 

A  mathematical  heat  transfer  model,  applicable  to  determining 
the  sensitivity  to  input  current  of  a  hot  wire  initiator,  is 
developed.  Simplifying  assumptions  are  introduced  which  make  the 
model  amenable  to  solution  using  a  medium  size  analog  computer. 

The  model  is  applied  to  the  MARK  45  MOD  0  impulse  cartridge.  Solu¬ 
tions  are  compared  to  empirical  sensitivity  data  for  this  cartridge 
conditioned  at  -30°F,  68  F  and  225° F. 
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INTRODUCTION 


A  hot  wire  initiator  (see  Figure  1}  is  a  single  discrete  unit 
containing  an  explosive,  propellant,  or  pyrotechnic  material 
commonly  referred  to  as  the  primary  mix  or  bead  material;  hereafter 
referred  to  as  the  bead.  The  bead  usually  surrounds  a  fine  wire 
known  as  the  bridgewire  which  spans  or  bridges  two  pins  as  shown  in 
Figure  1.  The  bead  is  ignited  by  ohmic  heating  of  the  bridgewire 
when  sufficient  electric  potential  is  applied  across  the  pins.  The 
ignited  bead  in  turn  ignites  a  secondary  charge  known  as  the  booster 
or  main  cha;  ge,  as  the  case  may  be . 

Of  interest  to  the  designer  of  a  hot  wire  initiator  is  the 
control  of  its  sensitivity.  Here  sensitivity  refers  to  the  relation 
between  the  firing  time  of  an  initiator  and  the  input  firing  stimulus 
such  as  current  or  potential.  The  problem  of  determining  the  sensi¬ 
tivity  of  a  newly  designed  hot  wire  initiator  was  treated,  in  the 
study  reported  here,  as  a  heat  transfer  phenomenon.  A  mathematical 
model  was  developed  whose  solutions  yield  the  sensitivities  of  hot 
wire  initiators.  The  model  was  then  applied  to  the  MARK  45  MOD  0 
Impulse  Cartridge.  In  o"der  to  obtain  values  for  two  inaccurately 
known  parameters  (the  thermal  conductivity  of  the  bead  and  the  heat 
transfer  film  coefficient  of  the  bridgewire-bead  interface),  they 
were  systematically  varied  until,  through  trial  and  error,  the  calcu¬ 
lated  sensitivity  matched  the  measured  sensitivity  at  two  values  of 
input  current.  In  this  way  the  model  was  calibrated  against  empirical 
firing  time  versus  firing  current  data  from  cartridges  conditioned  at 
room  temperature.  The  degree  of  fit  between  the  calculated  and 
empirical  data  for  other  values  of  input  current  and  at  different  con¬ 
ditioning  temperatures  was  used  as  a  measure  of  the  adequacy  of  the 
model. 

Finally,  certain  of  the  inaccurately  known  parameters  were 
arbitrarily  changed,  one  at  a  time,  to  determine  the  sensitivity  of 
the  solutions  to  possible  error  in  the  parameter  changed. 

All  solutions  were  obtained  on  an  Electronic  Associates  analog 
computer  consisting  of  an  eighty  amplifier  TR-10  four  cabinet  complex. 


MATHEMATICAL  MODEL 


X.  Equation  of  Conduction  of  Heat  - 

Reference  (a),  page  10,  gives  the  heat  conduction  equation  for 
an  isotropic,  homogeneous  solid  with  internal  heat  generation  as 


M  =aV2e  +  Ai ZjlZ&IL, 
bt  k 


(1) 
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where : 

a  =  iL.  is  the  thermal  diffusivity, 
pc 

A  is  the  rate  of  heat  produced  per  unit  time  per  unit, 
volume  at  the  point  (x,y, z), 

c  is  the  specific  heat, 

k  is  the  thermal  conductivity, 

p  is  the  density, 

6  is  the  instantaneous  temperature  at  the  point  x,y,z, 

^  is  the  Laplacian  operator. 

II.  Simplifying  Assumptions  - 

In  order  to  obtain  a  mathematical  model  within  reach  of  the  capa¬ 
bilities  of  the  computing  equipment  used,  the  following  simplifying 
assumptions  were  made; 

A.  Radial  temperature  distributions  in  the  bridgewire  were 
neglected.  Erosn  equation  (l),  page  204  of  reference  (a),  it  was 
concluded  that  the  difference  between  the  outer  surface  temperature 
and  the  center  temperature  of  the  bridgewire  of  the  MARK  45  MOD  0  car¬ 
tridge  (the  hot  wire  initiator  used  in  this  analysis)  would  never 
exceed  1.34° C.  This  difference  was  considered  negligible  since  it  is 
only  0 .5$  of  the  assumed  autoignition  temperature  of  the  bead. 

B.  Although  the  MARK  45  MOD  0  cartridge  is  a  four  pin  double 
bridge  hot  wire  initiator,  only  one  of  the  bridge  circuits  was  con¬ 
sidered  in  the  study.  Figure  2  is  a  sketch,  with  coordinate  system 
superimposed,  of  the  bridgewire  and  bead  of  the  MARK  45  MOD  0  car¬ 
tridge.  Only  one  of  the  bridgewlres  is  shown.  The  other  bridgewire 
is  parallel  to  the  one  shown  and  is  electrically  insulated  from  it. 

C.  Cylindrical  symmetry  was  assumed. 

D.  As  depicted  by  the  broken  lines  in  Figure  2,  an  arbitrary 
cyclindrical  boundry  was  introduced.  In  the  mathematical  model  this 
was  assumed  to  be  the  outer  boundary  of  the  bead  and  was  assumed  to 
remain  at  ambient  temperature.  This  assumption  was  based  on  the 
concept  that  the  temperature  of  the  relatively  large  masses  of  the 
pins  and  the  ceramic  plug  probably  does  not  differ  appreciably  from 
ambient . 
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E.  The  thermal  properties  (e.g.,  thermal,  conductivity; 
specific  heat,  and  diffusivity)  of  the  initiator  materials  were 
assumed  to  he  independent  of  temperature.  Approximately  median 
values  for  these  properties  were  chosen  consistent  with  the  tempera¬ 
ture  range  encountered . 

F.  The  exothermic  chemical  self  heating  of  the  bead  was 

assumed  negligible  until  the  bead  reached  a  critical  temperature  T  * 

a 

Ignition  was  assumed  to  occur  at  this  time  and  all  subsequent  calcula¬ 
tions  were  disregarded.  The  value  chosen  for  Ta  was  based  on  experi¬ 
mental  evidence  of  an  exothermic  reaction  in  the  bead  when  the 
initiator  was  subjected  to  various  oven  temperatures.  In  accord  with 
the  experimental  technique  for  determining  the  value  of  Ta,  it  is 
referred  to  as  the  autoignition  temperature. 

Ill .  Final  Equations  - 

Let  primed  symbols  refer  to  the  bridgewire  and  let  unprimed 
symbols  refer  to  the  bead .  Also  let  6  ( G ')  now  refer  to  the  tempera¬ 
ture  difference  between  the  bead  (bridgewire)  and  ambient.  As 
developed  in  sections  4.2  and  4.10  of  reference  (a);  equation  (1) 
applied  to  the  bridgewire  becomes 

dfl'  ,  I2  J(HaV)  P'h  fa*_a^  (2) 

TT - a  +  -77V  V'S-  /  '  / - 7T7-7  "sb 

ot  p(A)co  pAc 

where : 


A'  is  the  cross  sectional  area  of  the  bridgewire, 

a'  is  che  diffusivity  of  the  bridgewire, 

h  is  the  heal  transfer  film  coefficient  of  the  bridgevire- 
bead  interface , 

I  is  the  instantaneous  current  in  the  wire, 

J  is  the  mechanical  equivalent  of  heat, 
t  is  the  time , 

a'  is  the  temperature  coefficient  of  resistivity > 

Ss  is  the  surface  temperature  of  the  bead  at  the  bridgewire- 
bead  interface  . 
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Equation  (1)  applied  to  the  bead  simplifies  to 

||  =  s^e.  (3) 

where  a  is  the  diffusivity  of  the  bead. 

The  initial,  boundary,  and  terminal  conditions  ere 


9=q' -  0  '-t  t  =  0, 

(4) 

e-"-0/  =0  at  x  =  V  |,  r=rbcad. 

(5) 

d'(-x)  ^  ?'(x) , 

(6) 

e(-x)  =  5(x) , 

(7) 

-k  —  =  h(d'-0s)  at  r=rs, 
dr 

(8) 

9  <  Ta-Tc  , 

(9) 

where, 

L  is  the  length  of  the  bridgewire  and  the  bead, 

rbead  ou'^er  radius  of  the  bead, 

rs  is  the  inner  radius  of  the  bead  (equal  to  the  bridge- 
wire  radius), 

Tc  is  the  ambient  temperature . 

Equations  (2)  and  (3)  can  be  approximated  by  a  set  of  ordinary 
differential  equations  with  time  as  the  independent  variable.  This 
is  done  by  considering  9 '  at  discrete  values  of  x  and  by  considering 
9  at  discrete  values  of  x  and  r.  The  partial  derivatives  of  9 1  with 
respect  to  x  and  of  9  with  respect  to  x  and  r  are  then  represented  by 
appropriate  difference  expressions  (see  Appendix  A). 

Two  independent  investigations  were  made  to  determine  the  distri¬ 
bution  and  the  minimum  number  of  discrete  values  of  x  and  r  at  which 
the  temperatures  must  be  computed  to  obtain  reasonable  accuracy. 
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Computational  experiments  involving  up  to  6  axial  temperature  sta¬ 
tions  in  half  the  bridgewire  length  (equivalent  to  12  axial  sta¬ 
tions)  showed  that  3  equally  spaced  axial  stations  in  half  the 
bridgewire  and  bead  length  was  sufficient.  The  error  in  computed 
temperatures  introduced  by  using  only  3  axial  stations  was  about 
3$.  Similarly,  computational  experiments  involving  up  to  6  radial 
temperature  stations  showed  that  4  unequally  spaced  radial  stations 
in  the  bead  were  sufficient.  The  error  in  computed  temperatures 
introduced  by  using  only  4  radial  stations  was  less  than  10$.  This 
was  provided  the  radial  temperature  station  nearest  the  bridgewire- 
bead  interface  was  chosen  sufficiently  close  to  the  interface.  It 
was  determined  that  "sufficiently  close"  meant  that  the  distance 
from  the  bridgewire-bead  interface  to  the  first  radial  temperature 
station  should  be  less  than  1$  of  the  distance  from  the  bridgewire- 
bead  interface  to  the  bead's  outer  radius. 


In  accordance  with  the  conclusions  reached,  temperature  sta¬ 
tions  were  selected  within  the  bridgewire  and  bead  of  the  MARK  45 
MOD  0  cartridge  as  indicated  in  Figure  3.  Ea-h  station  is 
assigned  the  generic  indices  n,m  where  n  is  the  axial  position 
index  and  m  is  the  radial  position  index,  n  =o  refers  to  any  point 
in  a  plane  perpendicular  to  the  bridgewire  and  passing  through  its 
center.  m=o  refers  to  any  point  on  the  surface  of  a  circular 
cylinder  whose  axis  coincides  with  the  axis  of  the  bridgewire  and 
whose  radius  is  equal  to  the  radial  distance  from  the  axis  of  the 
bridgewire  to  the  row  of  stations  within  the  bead  nearest  ihe 
bridgewire . 

For  the  bridgewire 


n2  „ 
o  a 


S?" 


(1C) 


From  equation  (A  13)  in  Appendix  A,  equation  (ID)  can  be  approxi¬ 
mated  (with  p=lj  by 


20'  ♦ 


(11) 
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Substituting  from  equation  (11)  into  equation  (2),  combining  terms, 
and  imposing  the  initial,  boundary,  and  symmetry  conditions  of 
equations  (4),  (5),  (6)  and  (7),  the  final  equations  for  the  bridge- 

wire  are 


In  order  to  calculate  GU}  5,  a  parabola  lying  in  the  r -6  plane 
is  fitted  through  the  three  points  (-Po>0n..i)>  (O,0njO)j  and 
(?Po#0n  x)  818  shown  in  Figure  4.  The  equation  of  this  parabola 
with  origin  of  the  coordinate  system  arbitral ily  chosen  at  radial 


station  r= 


9  A(r’r0'2  +  B(r-r0)  +  C* 


9n,s  =Al"  7T 


Tte  coefficients  A,  B,  and  C  are  derived  in  Appendix  A.  Substituting 
the  values  for  A,  B,  and  C  from  equations  (AS),  (A9),  and  (A10)  into 
equation  (14)  results  in  the  final  expressions  for  the  bead  inner 
surface  temperatures. 

0_  =  -  v1-.  ■  0_  ,  +  0_  .  +  6_  ,  1 


4 7+1 )  n>  i  4  7  n>  0  4~(  7+1)  n, -1 


=  0,1,2. 
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The  final  form  of  the  boundary  condition,  equation  (8),  can 
now  be  given.  Substituting  into  equation  (8)  the  finite  difference 

approximation  for  derived  in  Appendix  A  (see  equation  (All)}, 

a..d  the  expressions  for  the  bead  surface  temperature  (equation  (15)), 
and  combining  terms  results  in 


9n,  -1 


1  h 
47(V+i)  k 


- * -  01  +  ■  ./ V  /LTtL,  -  6n  i  + 

i_  +  2z±L  u  +  2y+ i  h 


h  ,  2-y+l  h  11,1  1.  2  7+1  h 

.k  4(7+1)  k /  \pn  4(7+1)  k 


27+I  h\ 

24li"h  rn'c 


(16) 


n  =  0,1,2. 


For  the  bead 


y 20  =  +  iL£. 

dr2  r  dr  dx2 


From  equation  (A14)  in  Appendix  A,  equation  (17)  can  be  approximated 
(with  p=l)  by 


V20  x=3Cv 


2r  +p 
m  m 


r^n  \r(7+l)P2rm 


*n,  m+1 


( 2r  +(1-7)  p  0 
I  m  ''  -^m  A  2  ,  n 

- __ -  +  ~2  yn,  m 


+  U2  y 
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Substituting  from  equation  (IS)  into  equation  (3)  and  imposing  the 
initial,  boundary,  and  symmetry  conditions  of  equations  (4),  (5), 
(6),  and  (7);  the  final  equations  for  the  bead  are 


Equations  (12),  (15),  (16),  and  (19)  are  the  final  equations  for 
the  bridgevire,  bead  inner  surface  temperature,  boundary  condition, 
and  bead,  respectively. 
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RESULTS 


Enpirical  data  of  firing  time  versus  constant  DC  current, 
for  the  MARK  45  MCD  0  cartridge  at  three  conditioning  tempera¬ 
tures,  is  presented  in  Figure  5.  The  vertical  bars  indicate  the 
spread  in  the  data.  A  smooth  curve  (not  shown)  was  visually 
fitted  to  the  room  temperature  data,  and  from  this  curve,  the 
sensitivity  at  4  and  10  amps  was  determined.  Several  values  for 
the  thermal  conductivity  k  of  the  bead  were  chosen  which  were 
expected  to  bracket  the  true,  but  unknown,  value.  Using  currents 
of  4  and  10  amps  in  the  computer  program,  selected  values  of  heat 
transfer  .film  coefficients  h  for  the  bridgewire-bead  interface 
were  determined  by  trial  and  error.  These  values  were  selected 
such  that  the  computed  sensitivity  matched  the  empirical  sensitiv¬ 
ity  for  each  value  of  k  chosen.  By  graphing  the  corresponding 
values  of  h  vs  k  for  4  amps  and  for  10  amps  respectively,  it  was 
determined  that  a  single  value  of  h  and  of  k  would  satisfy  the 
empirical  sensitivity  requirement  at  both  4  and  10  amps.  This 
result  is  shown  in  Figure  6.  This  indirect  method  of  measuring  h 
and  k  is  dependent  on  the  adequacy  of  the  mathematical  model  and 
the  accuracy  of  the  other  parameters.  Due  to  the  large  number  of 
assumptions  required  in  this  study,  no  attempt  was  made  to  estimate 
the  accuracy  of  the  values  obtained  for  h  and  k. 

Table  1  lists  the  input  parameters  used  in  these  and  subsequent 
calculations.  It  also  lists  the  composition  of  the  initiator 
materials  of  interest  here  and  the  sources  for  information  given. 

TABLE  1  INPUT  PARAMETERS  AND  MATERIAL  COMPOSITION 


Parameter  Symbol 

Bridgewire : 

Length  L,  6q 

Radius  r 

Cross  sectional  area  A ' 

Perimeter  P ' 

Electric  conductivity  o' 

Thermal  conductivity  k ' 

Density  p ' 

Thermal  capacity  c ' 

Temperature  coefficient  a! 

of  resistivity 


Value  Source* 


0.2847  cm  1 
0.9010  x  10 -2  cm  1 
2.554  x  10 -4  cm2  1 
5.661  x  IQ-2  cm  1 
1.00  x  104  mhos/cm  2 
0.050  cal/sec-cm-°C  3 
8 . 50  gm/ cm3  4 
0.126  cal/gm-° C  4 
1.8  x  10"4  (°C)_1  2 
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TABLE  1  (Continued) 


Parameter 

Symbol 

Value 

Source* 

Bead: 

Length 

L,  6q. 

0.2847  cm 

1 

Inner  radius 

rs 

0.9010  x  10 ~2  cm 

1 

Outer  radius 

rbead*  r4 

7.875  x  10 ~z  cm 

1 

Thermal  conductivity 

k 

5.05  x  10 ~4  cal/sec-cm-°C 

5 

Density 

P 

2.52  gm/cm3 

4 

Thermal  capacity 

c 

0.2461  cal/gm-°C 

4 

Autoignition 

270°  C 

6 

temperature 

Brldgevire-bead  Interface: 

Heat  transfer 

h 

0.116  cal/sec-cm2-°C 

5 

film  coefficient 

Bridgewire  composition: 

Nickel 

Ni 

80$ 

1 

Chromium 

Cr 

20$ 

1 

Bead  composition: 

Zirconium 

Zr 

7 .00$ 

1 

Ammonium  dichromate 

(NH4)2Cr207 

14.08$ 

1 

Ammonium  perchlorate 

kh4cio4 

44.92$ 

1 

Barium  chromate 

BaCr04 

34.00$ 

1 
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TABLE  1  (Continued) 

♦Source  key 

1.  Either  obtained  directly  or  derived  from  information 
given  in  reference  (d). 

2.  Obtained  directly  from  reference  (e). 

3.  Derived  from  information  given  in  reference  (f). 

4.  Derived  from  information  given  in  reference  (e).  The 
derivation  was  based  on  values  given  for  each  component  of  the 
composition  if  it  was  available.  If  it  was  not  available,  the 
values  for  chemically'  similar  compounds  were  used. 

5.  Derived  by  graphically  fitting  computed  results  to 
empirical  data  at  two  input  currents. 

6.  Based  on  experimental  evidence  of  an  exothermic 
reaction  in  the  bead  when  the  initiator  was  subjected  to  various 
oven  temperatures . 


Using  the  values  of  k  =  5.05  x  10 “4  cal/cm-° C  and  h  =  0.116 
cal/cm2-°C  (the  point  of  intersection  in  Figure  6),  the  bead  tempera¬ 
ture  0O , s  of  the  MARK  45  mop  0  cartridge  was  calculated  for  values 

of  input  current  ranging  from  3.0  to  20  .0  amps.  The  o,s  point  was 
selected  because  it,  being  in  the  center  and  at  the  bead's  inner  sur¬ 
face  (nearest  the  bridgewire),  is  always  the  hottest  point  in  the 
bead.  The  initiator  sensitivity  for  a  given  current  and  condition¬ 
ing  temperature,  was  then  determined  by  observing  the  time  required 
for  Gc j s  to  reach  a  value  equal  to  the  difference  between  the  assumed 
autoignition  temperature  of  ?70°G  and  the  conditioning  temperature. 
Conditioning  temperatures  of  -80,  68,  and  225°F  were  investigated. 
These  correspond  to  -62.22,  20.00,  and  107. 2° C,  respectively.  The 
computed  sensitivities  at  the  three  conditioning  temperatures  are 
presented  by  the  three  curves  in  Figure  5. 

Figure  7  is  a  plot  of  90}  s  versus  time  for  various  currents. 

The  intersections  of  these  curves  with  the  horizontal,  lines  gives 
the  times  at  which  ignition  is  assumed  to  occur.  The  sensitivity 
curves  of  Figure  5  are  derived  from  this  infoimation. 


In  Figure  7,  the  curve  of  0Oj  5  versus  time,  for  1=3  amps, 

indicates  that  a  steady  state  temperature  below  the  autoignition 
temperature^ would  be  reached  at  condit?  ming  temperatures  below 
-10° C  (14.0  F).  Hence,  as  an  example,  steady  state  temperatures 
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were  computed  for  the  bridgewire  and  bead  for  1=3  amps.  Figure  8 
is  a  plot  of  these  steady  state  temperatures  versus  radius  r. 

Figure  9  is  a  plot  of  the  same  steady  state  temperatures  versus 
axial  position  x. 

Table  1  shows  tliat  values  for  h,  k,  p,  and  c  (not  being 
available  in  the  literature)  were  estimated  from  related  data.  The 
geometry  of  the  outer  boundary  was  arbitrary  and  the  autoignition 
temperature  Ta  was  based  on  a  gross  experiment.  To  determine  the 

sensitivity  of  the  solutions  to  each  of  these  parameters,  the  sen¬ 
sitivity  of  the  MARK  45  MOD  0  cartridges  was  recalculated  with  each 
of  the  above  parameters  changed  one  at  a  time.  Since  p  and  c  appear 
only  as  a  product,  only  the  product  pc  was  changed.  The  input 
parameters  Ta,  h,  and  k  were  each  increased  20$.  The  product  pc  was 

decreased  20$.  The  bead  radius  and  the  bridgewire  and  bead  length 
were  each  increased  50$.  Each  of  these  changes  is  reflected  in  the 
sensitivity  of  the  initiator.  In  Figures  10  through  15,  the  new 
sensitivity  is  compared  to  that  originally  calculated  and  shown  in 
Figure  5.  The  comparisons  are  made  for  a  conditioning  temperature 
of  68°F. 


CONCLUDING  REMARKS 

A  mathematical  heat  transfer  model  has  been  developed  which  is 
amenable  to  solution  using  a  medium  size  analog  computer.  The  model 
is  applicable  to  determining  the  sensitivity  of  a  conventional  hot 
wire  initiator.  While  the  model  has  been  applied  only  to  the  MARK  45 
MDD  0  impulse  cartridge,  its  ability  to  faithfully  reproduce  the  room 
temperature  sensitivity  of  this  cartridge  over  a  range  of  firing 
times  from  15  to  1000  milliseconds  demonstrates  the  adequacy  of  the 
model  to  account  for  the  important  parameters  associated ’with  a  hot 
wire  initiator  of  this  basic  design.  Additionally,  the  model  predicts 
with  reasonable  accuracy  the  sensitivity  of  the  MARK  45  MOD  0  impulse 
cartridge  conditioned  at  22 5° F  and  at  -80°  F. 

At  present,  fundamental  data  required  as  input  to  the  model  is 
only  partially  available.  However,  the  usefulnes  .  of  the  model  is 
not  lost  since,  as  was  done  in  this  study,  the  model  can  be  calibrated 
to  reproduce  the  available  sensitivity  data.  The  sensitivity  for  other 
conditions  can  then  be  predicted.  For  example,  the  model  is  not 
restricted  to  DC  input  currents.  However,  if  the  sensitivity  to 
DC  input  currents  is  known  then  the  model  can  be  calibrated  to  repro¬ 
duce  this  data.  Subsequently,  the  model  could  serve  as  a  powerful  t col 
to  predict  the  sensitivity  to  any  time  varying  input  current,  e.g. 
sinusoidal.  Or,  it  could  be  used  to  predict  (or  investigate)  the 
effects  on  sensitivity  of  changing  one  or  more  design  parameters  such 
as  the  bridgewire  material  or  bridgewire  diameter. 
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The  usefulness  of  the  model  as  a  starting  point  for  the 
design  of  nev  hot  wire  initiators  is  yet  to  be  demonstrated.  How¬ 
ever,  the  work  reported  here  indicates  the  model  is  adequate  to 
predict  representative  firing  times  versus  input  stimuli  of  conven¬ 
tional  hot  wire  initiators. 
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APPROXIMATION  OF  FIRST  AND  SECOND  PARTIAL  DERIVATIVES  AND  THE  LAPIACIAN 
BY  FINITE  DIFFERENCES  -  UNEQUAL  SECMENTS  AND  CYLINDRICAL  SYMMETRY 


The  Laplacian  in  cylindrical  coordinates  for  a  system  having 
cylindrical  symmetry  is 

v2@  =i££  +1^9  +  (Al) 

dr2  r  dr  dx2 

d2g 

To  evaluate  — and  a  parabola  lying  in  the  r-0  plane  with  x=x_  is 
dr  dr  n 

fitted  through  the  three  points  (-pffi, 0n^  n-1),  (0,0^  m),  and  (xPm>0n,m+1) 

as  shown  in  Figure  16.  The  equation  of  this  parabola  with  origin  of 

the  coordinate  system  arbitral  ily  chosen  at  radial  station  r=*r  is 

m 


Also, 


x=oc_ 


=  A( r-r  )2  +  B(r-r  )  +  C, 


(A2) 


d0 

3? 


=  B, 


x=aCj 

r=r, 


(A3) 


£0 

dr 


=  2A , 


X=0Cj 

r=r 


ra 


(A4) 


0 


n,  ra-1  =A(-P„)2  4  +  C’ 


°„,m  =A(°)2  +  ®(0)  +  C- 


(AS) 

(A6) 


and 


n,  n+i 


-Mzp.r 


B(7Pm) 


+  C. 


(A7) 
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Simultaneous  solution  of  equations  (.AS),  (AS),  and  (A7)  for  A,  B, 


and  C 


A  = 


P®  7(7+1) 


9n,a«  -  (nl)  9n,n  •: 


(AS) 


B  = 


pra  7(7+!)  /n>m+1  '  (1“72)  en,m  -  7sen>a_1 


(A9) 


C  =  0n,m' 


^SSlne4UEtl0nS  (A8)  (A9)  *ith  *«“«<»•  («)  »d  (A4) 


(A10) 


M 

dr 


x=xr 

r=r 


7(7+1)  pm  L  n^+1 


(I-72)  9  -  yzQ 

r,m  n,a~3 


m 


(All) 


£0 

dr2 


x=xr 

r=r 


7(7+1)  p2  l8"-”+l  '  (l47)  9n,m  +  ^S„,a-1 


m 


m 


(AZ2) 


To  evaluate  -f  a  parabola  lying  in  the  x,  e  plane  with  r*r  is  fitted 
through  the  three  points  (~<ln>0  j  ),  (0,0  )  and  (*  ,e" 

shown  in  Figure  16.  Analogous  to  ihe  approbation  oh^edbr 


£2q 

;}  a  similar  approximation  for  v  e  is 

dx2 


dr2' 


dfe 

dx2 


x=oc 

r 

r=r_ 


0(1+6) 


‘n 


?n+l,m  “  (1*6)  0n,m  +  0  0n-l,; 


(A13) 
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Combining  equations  (All),  (A12),  and  (A13)  with  equation  (Al)  results 
in 
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